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ABSTRACT 

We present the results of a near-infrared monitoring program of the Anoma- 
lous X-ray Pulsar IE 2259-1-586, performed at the Gemini Observatory. This 
program began three days after the pulsar's 2002 June outburst, and spans ~1.5 
years. We find that after an initial increase associated with the outburst, the 
near-infrared flux decreased continually and reached the pre-burst quiescent level 
after about one year. We compare both the near-infrared flux enhancement and 
its decay to those of the X-ray afterglow, and find them to be remarkably consis- 
tent. Fitting simple power laws to the RXTE pulsed flux and near- infrared data 
for t > 1 day post-burst, we find the following decay indices: a = —0.21 ± 0.01 
(X-ray), a = —0.21 ± 0.02 (near- infrared), where flux is a function of time such 
that F (xf^. This suggests that the enhanced infrared and X-ray fluxes have a 
physical hnk post-outburst, most likely from the neutron-star magnetosphere. 

Subject headings: pulsars: individual (IE 2259-1-586) — pulsars: general — in- 
frared: stars — stars: neutron — stars: magnetic fields 



1. INTRODUCTION 

The origin of the optical and infrared (IR) emission in Anomalous X-ray Pulsars (AXPs) 
is currently uncertain. To date, possible IR counterparts have been identified for 5 (of 6) 
known AXPs [IE 2259-^586: HuUeman et al. (2001); Kaspi et al. (2003), IE 1048.1-5937: 
Wang & Chakrabarty (2002); Israel et al. (2002), IRXS J170849-400910: Israel et al. (2003), 
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4U 0142+61: HuUeman et al. (2004) and XTE J1810-197: Israel et al. (2004)]. In one case, 
4U 0142+61, optical emission has been detected (Hulleman et al. 2000a) and was seen to 
pulse at the same period as the X-ray pulsar (Kern & Martin 2002). All of these sources 
show excess optical/lR emission when compared to the prediction of a simple blackbody 
model extrapolated from X-ray energies (assuming the 2-component model consisting of a 
power-law plus blackbody component for the X-ray emission). The extrapolation of the 
power-law component greatly overpredicts the optical/IR flux, however. 

On 2002 June 18, IE 2259+586 exhibited an outburst that included, apart from the 
>80 bright, short-lived X-ray bursts, a large and long-lived X-ray flux enhancement with 
subsequent decay (Kaspi et al. 2003; Woods et al. 2004; Gavriil et al. 2004), as well as a 
signiflcant near-IR flux enhancement, demonstrated with Gemini-North Target of Opportu- 
nity observations made 3 and 10 days after the outburst (Kaspi et al. 2003). This was the 
flrst demonstration of IR variability in an AXP, and subsequently the flrst example of an 
associations between X-ray and IR activity. Since then, IR variability has also been reported 
in IE 1048.1-5937 (Wang & Chakrabarty 2002; Israel et al. 2002; Durant et al. 2004) and 
4U 0142+61 (Hulleman et al. 2004) without evidence for a correlation with X-ray outbursts, 
though sparsity in X-ray observations do not preclude this. Recently, Rea et al. (2004) re- 
ported that a second near-IR observation of the proposed counterpart to XTE J1810— 197 
showed variabhty possibly hnked to X-ray flux decay. 

In order to verify that the enhancement in IR flux seen post-outburst in IE 2259+586 
was genuinely related to the outburst, as well as to characterize its decay, we monitored the 
source using the Gemini North telescope. Here we report on this program, demonstrating 
conclusively that the IR enhancement reported by Kaspi et al. (2003) was associated with 
the outburst. We also flnd that the post-outburst IR and X-ray radiation properties are 
correlated. In §4, we compare our results with expectations from different models and other 
AXPs. 



2. OBSERVATIONS 

Images in the near-IR Kg band (A = 2.15 /xm, AA = 0.31 /xm) were obtained with 
the Near-InfraRed Imager (NIRI; f/6 camera; Aladdin InSb detector array; 1024 x 1024 
pixels; 0'.'1171 pixel size) at the 8-m Gemini North Observatory. For a description of the 
observing parameters and conditions, see Table 1. The detector array was read several times 
in order to reduce read-out noise. This was done both before and after each exposure, and 
the difference was recorded in the data files. Each target frame consisted of 4 coadded 
exposures of 15-s integrations. Two earlier observations of IE 2259+586, which took place 
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3 and 10 days after the X-ray burst in June 2002 as part of a Target of Opporunity (ToO) 
program, were described in Kaspi et al. (2003) and are also included in this analysis. All 
data were reduced using the Gemini IRAF package for NIRI data. Each frame was divided 
by a normalized flat fleld constructed from the lamp flat frames obtained with the Gemini 
standard calibration unit. The sky background image was derived from the data frames 
themselves with the objects masked out, taking advantage of the 9-point dither pattern 
applied during the observation: this was subtracted from all data frames. Finally, all data 
from a single night were coadded into one image. 

Photometry was performed using standard procedures within the DAOPHOT package 
(Stetson 1987) as implemented in IRAF. To calibrate the instrumental magnitudes found by 
DAOPHOT, we tied our measurements directly to stars in the IE 2259+586 field (HuUeman 
et al. 2001). The eight nearest neighbours that were bright, isolated and not varying (stars 
A, B, B', D, F, G, K and N under the numbering system of HuUeman et al. 2001) were used 
to measure the mean offset between the instrumental and published Ks band magnitudes. To 
verify non- variability, we chose only stars that fell within a standard devation of cr < 0.03 mag 
from the weighted mean. The advantage of this procedure was that it gave much more precise 
relative fluxes than would have been possible by using the single standard star observed 
during each night. 

Measured Ks band magnitudes of IE 2259+586 are listed in Table 1. We analytically 
estimated magnitude errors from the standard deviation of the sky background, under the 
assumption that an aperture with a radius equal to the PSF FWHM contains 70% of a star's 
flux: this was added in quadrature to the uncertainty in the photometric tie to produce flnal 
uncertainties in Table 1. Our careful re-analysis of the ToO data (the first two points) gave 
results consistent with Kaspi et al. (2003) well within uncertainties: the differences in Kg 
magnitude between the two analyses were 0.05 ± 0.17 (Jun. 21) and 0.18 ± 0.25 (Jun. 28). 
The third data point was observed with the Canada France Hawaii Telescope (CFHT) by 
Israel et al. (2004): we have assumed K' = Kg in our magnitude to flux conversion. From 
Table 1, one sees that ~400 days post-burst, the source appears to have returned to its 
pre-burst brightness of Kg = 21.7 ± 0.2 mag (HuUeman et al. 2001). 



3. RESULTS 

The post-outburst evolution of the X-ray pulsed flux of IE 2259+586 is described by 
Woods et al. (2004) in terms of a model with two power laws in time, where F cct", with F 
the unabsorbed 2—10 keV pulsed flux, t the time since the glitch epoch tg = 52443.13 MJD, 
and a the power-law index. Immediately after the outburst (<1 day) the decay appears to 
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foUow a much steeper power-law index than it does over the following year. We compare the 
near-IR flux enhancement and decay to those of the second, slower X-ray segment consisting 
of flux from >1 day post-burst, for which Woods et al. (2004) find a temporal decay index 
of a = -0.22 ±0.01. 

To the near-IR data, we first apply an extinction correction of Ak, — 0.56 ±0.01, which 
is inferred from Ay = Nh/{1.79 x 10^^ cm~^) = 5.2 mag (Predehl & Schmitt 1995) where 
Nh = 9.3 ± 0.3 X 10^^ cm-2 (Patel et al. 2001). Magnitudes are then converted to uF^ (see 
Table 1 and Figure 1). 

We fit a simple power-law function, F = A;(t/100)°, to the X-ray and IR flux using 
a numerical fitting program that directly searches over parameter space, where /c is a 
constant with dimensions erg/s/cm^, t is time in days since the glitch, and the choice of the 
factor of 100 roughly minimizes the covariance between k and a. Note also that F represents 
both X-ray flux and near-IR pF^, depending on the case, in erg/s/cm^. Table 2 contains 
our best-fit parameter values and la uncertainties. Fitting the RXTE pulsed-flux data only, 
we confirm the index a reported by Woods et al. (2004)^. Comparing the long-term decay 
of near-IR and X-ray afterglow, we find the simple power-law indices remarkably consistent: 
a = -0.21 ± 0.01 (X-ray) and -0.21 ± 0.02 (IR). 

Interestingly, the X-ray and IR enhancements also appear to be offset from their respec- 
tive quiescent levels at t = 1 by nearly the same amount and, as a result, decay back to Fq 
on similar time scales. To quantify this, we perform a fit to a second function consisting 
of a power-law with an excess offset F — Fq{l ± /(i/io)°), where to — 3.5 days is the time 
since the glitch of the first IR observation, / = [Fq — Fq)/Fq is the flux excess, and Fq is 
the flux enhancement at to- Data from t < are included to determine Fg. If the offsets are 
correlated, then we would expect the best-fit / from X-ray and Kg data to be consistent: 
this is in fact what we find in the latter part of Table 2. We note that larger errors and a 
small number of near-IR data points result in x^/ v < 1 and large uncertainties on / and 
a when all three parameters are fit to the near-IR data. To confirm that the shapes of the 
decay curves of the X-ray and Kg data are consistent statistically, we re-fit the Kg data to 
the excess model with a and / fixed at the X-ray excess best-fit values. The effect of holding 
two parameters constant is a reasonable increase in /v to 5.7/5 ~ 1.1, showing that the 
two shapes are indeed consistent with each other at the la level. In Figure 1, we plot the 
best- fit power laws as modelled on the X-ray and Kg data (dashed lines); overplotted is the 



^Thc small discrepancy can be explained by noting that Woods et al. (2004) perform a fit on log- 
log data to a linear function, which neglects to account for asymmetric uncertainties, unlike our numerical 
method. 
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power-law plus excess model where a — —0.44 and / = 2.14 (dot-dashed lines), and the 
corresponding Fg best- fit value (dotted lines), for comparison. 

Thus, we find that the >1 day IR and X-ray initial enhancements and subsequent decays 
after the 2002 outburst are correlated. 



4. DISCUSSION 

In this section, we compare our results with other AXPs and discuss them in the context 
of various models. 

AXP optical/IR emission has been argued as originating from a fossil disk around the 
neutron star (e.g. Chatterjee et al. 2000; Eksi & Alpar 2003). The exceptionally high X-ray to 
optical/IR flux ratio seen in AXP IE 2259+586 rendered this model problematic (HuUeman 
et al. 2000b; HuUeman et al. 2001). Moreover, the SGR-hke bursting phenomena, observed in 
the 2002 outburst of IE 2259+586 (Kaspi et al. 2003; Gavriil et al. 2004; Woods et al. 2004), 
and from IE 1048.1-5937 in 2001 (Gavriil et al. 2002) and very recently in 2004 (Kaspi et al. 
2004), as well as the high optical pulsed fraction seen in AXP 4U 0142+61 (Kern & Martin 
2002), simply cannot be explained by such disks. Ertan & Cheng (2004) argue that the high 
optical pulsed fraction can be reproduced in a disk-dynamo model, although whether this 
could produce the observed bursts is unclear. 

"Hybrid" fallback disk models, in which the disk surrounds a magnetar, have recently 
been proposed to attempt to explain all AXP properties (Eksi & Alpar 2003). In this 
case, the quiescent pulsed X-ray emission arises from accretion from the disk while optical 
pulsations and bursts are magnetar magnetospheric emission. Ertan & Alpar (2003) argue, 
in the context of such a hybrid model, that enhanced X-ray emission, following an SGR-like 
flare, is released from the inner disk which has been pushed back by the burst itself. For 
IE 2259+586, however, no such SGR-like flare was detected prior to the observed enhanced 
X-ray emission, with upper limit 3 orders of magnitude below the total observed energy 
release (Woods et al. 2004). Furthermore, one AXP shows uncorrelated torque and X-ray flux 
variations, contrary to the predictions of any fallback disk model (Gavriil & Kaspi 2004). It 
is true that one natural prediction of fallback disk models is some form of correlation between 
the IR and X-ray emission (see Rea et al. 2004, and references therein). Nevertheless, in the 
absence of solutions to the problematic aspects of these models as listed above, we do not find 
the observed IR/X-ray correlation to render fallback disk models particularly compelling. 

By contrast, the magnetar model accounts very well for the bulk of AXP properties, 
especially bursts (Thompson & Duncan 1996), and qualitatively can explain optical/IR 
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properties as well. In the magnetar model, thermal surface emission is ruled out as the 
energy source for AXP optical/IR emission because of the impossibly high implied bright- 
ness temperature; hence, optical/IR emission must originate in the stellar magnetosphere, 
regardless of what powers it. Recently, Eichler et al. (2002) and Ozel (2004) argued that 
electron/positron radiation in the magnetosphere of a magnetar, produced in analogy with 
that in rotation-powered pulsars, could explain the observed optical/IR properties of AXPs. 
The strong correlation between IR and X-ray flux decay that we have observed in the after- 
glow of the 2002 outburst of IE 2259-1-586 argues for a physical link between the origins of 
both types of radiation. The X-ray emission is far too luminous {Lx ~ 10^^ in 1—10 keV at a 
distance ~4 kpc, Mereghetti et al. 2002) to be rotation-powered; it follows that the IR emis- 
sion is likely not cither. Though Ozel (2004) argued that the IR emission from an energetics 
standpoint could be rotation-powered, the large implied efficiency [vF^y / vFy^j-ot ~ 0.6, from 
Figure 1 of Ozel 2004) of conversion from spin-down flux into IR emission in 4U 0142+61, if 
the latter is rotation-powered, also argues against (though does not disprove) this hypothesis. 

The post-outburst correlation does clarify, however, the origin of the X-ray afterglow 
following the 2002 outburst, because of the following reasoning. Woods et al. (2004) identifled 
two possible mechanisms to produce the X-ray afterglow. The flrst was a genuine afterglow, 
i.e. thermal emission from the surface, a result of an impulsive heat injection to the crust 
from the magnetosphere. Such a thermal afterglow mechanism has been invoked in SGRs, in 
which the impulse, namely a bright soft gamma-ray flare, was clearly observed (Hurley et al. 
1999). For IE 2259-^586, however, no such flare was seen (Woods et al. 2004). An alternative 
to this thermal afterglow model is that the enhanced X-rays are a result of a twisting of the 
magnetosphcric field, perhaps as a result of the twisting of its footpoints following a signficant 
surface restructuring. Such an event is consistent with the coincidental rotational glitch that 
was observed (Kaspi et al. 2003; Woods et al. 2004), since the latter clearly implies a major 
disturbance in the crustal superfiuid. Such a twisting could naturally result in enhanced 
X-rays (Thompson et al. 2002), with subsequent field relaxation accounting for the decay. 
Given that the IR enhancement cannot be from surface thermal emission, the correlation 
with the X-ray decay strongly favors the twisting model, as it is difficult to understand how 
surface thermal X-ray emission could be so closely correlated with magnetosphcric radiation. 
Moreover, the IR enhancement being a result of a decaying magnetosphcric disturbance is 
consistent with the picture suggested by Eichler et al. (2002) and Ozel (2004) in which the 
IR emission is radiation from magnetosphcric pairs. It would be interesting to see, in future 
outbursts, whether the IR emission is pulsed and/or polarized, and if so, whether the pulse 
morphology is similar, and changes in concert with any X-ray pulse morphological changes, 
as this would strongly support this scenario. 

IR variability over long time periods has been seen in three other AXPs. In the case 
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of IE 1048.1-5937, X,-band variability has been detected (Wang & Chakrabarty 2002; 
Israel et al. 2002). However, no variation in X-ray flux was seen between the different 
epochs. Moreover, a more recent observation (Durant et al. 2004) found that the IR flux 
was consistent with the fainter of the two previous measurements, even though the X-ray 
flux was signiflcantly larger (Gavriil & Kaspi 2004). Despite showing no evidence of X-ray 
activity (Gavriil & Kaspi 2002, and unpublished work), 4U 0142-1-61 also appears to vary in 
Kg (HuUeman et al. 2004). This suggests that the physical mechanism responsible for the 
IR emission is distinct from that responsible for the quiescent X-rays, even during the broad 
X-ray flaring reported by Gavriil & Kaspi (2004). This is in contrast to the post-outburst 
behavior we have seen in IE 2259-1-586. That AXP IR and X-ray emission is generally 
correlated following an AXP outburst is further supported by the recent report of IR decay 
in data from the transient AXP candidate XTE J1810-197 (Rea et al. 2004) following its 
X-ray brightening and fading in 2004 (Ibrahim et al. 2004). 

This work was based on observations obtained at the Gemini Observatory (Program IDs 
GN-2002A-DD-6, GN-2003A-Q-71, GN-2003B-Q-22), which is operated by the Association 
of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF 
on behalf of the Gemini partnership: the National Science Foundation (United States), the 
Particle Physics and Astronomy Research Council (United Kingdom), the National Research 
Council (Canada), CONICYT (Chile), the Austrahan Research Council (Austraha), CNPq 
(Brazil) and CONICET (Argentina). It was also supported by NSERC Discovery Grant 
228738-03, NSERC Steacie Supplement 268264-03, a Canada Foundation for Innovation 
New Opportunities Grant, and FQRNT Team and Centre Grants. V. M. K. is a Canada 
Research Chair and Steacie Fellow. 
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Table 1. Near-IR observing parameters and measured results 



Date 




Instrument 


Exposure 
(min) 


Seeing 


Band 


Absorbed 
magnitude 


Unabsorbed vFy'^ 
(10-15 erg/s/cm^) 


2002 Jun. 


21 


Gemini/NIRI 


19 


O'.'T 


Ks 


20.41(7) 


10.7 ±0.7 


2002 Jun. 


28 


Gemini/NIRI 


12 


O'.'S 


Ks 


20.96(14) 


6.4 ±0.8 


2002 Aug. 


18 


CFHT/AOB/KIR 


122 


0'.'2 


K' 


21.31(24) 


4.6 ± 1.0 


2003 Aug. 


11 


Gemini/NIRI 


58 


O'.'S 


Ks 


21.66(11) 


3.4 ±0.3 


2003 Nov. 


5 


Gemini/NIRI 


49 


0'.'3 


Ks 


21.54(5) 


3.8 ±0.2 



^vF^ = 9.28 X 10-'^ erg/s/cm2 for absorbed Ks = Q (Chapter 7 of Cox 2000). 



Table 2. Results of power-law and excess function fitting 



Function 


Fq (erg/s/cm^) 


k 


(erg/s/cm 




/ 


a 






X-ray PL 




(2.35 


±0.02) X 


10-11 




-0.21 ±0.01 


37.P 


1.0 


IR PL 




(5.02 


± 0.25) X 


10-15 




-0.21 ±0.02 


7.1*^ 


1.0 


X-ray excess 


(1.56 ±0.03) X 10-" 








2.14 ±0.14 


_n 44+0-02 


97.3^ 


1.0 


IR excess 


(3.49t°:22) X 10-15 








Z.UO_Q 25 




1.4 


0.5 



^Errors scaled to infer uncertainties on the parameters, values reflect those before rescaling. 
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Fig. 1. — Unabsorbed X-ray flux and near-IR vFy decay of IE 2259+586 as a function of 
time. RXTE pulsed flux data (from Woods et al. 2004, Figure 13) are represented by circles 
and refer to the left axis, Gemini and CFHT data (Table 1) are represented by squares and 
refer to the right axis. Best-fit power laws to the X-ray and near-IR data are shown in dashed 
lines. The power-law plus excess model with a and / fixed at best-fit X-ray values is shown 
in dot-dashed lines; the dotted lines denotes the fiux levels during quiescence as determined 
by the excess fit. 



